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ABSTRACT 

We investigate mechanisms of energy extraction from a rotating black hole in 
terms of negative energy- at- infinity. In addition to the Penrose process through 
particle fission, the Blandford-Znajek mechanism by magnetic tension, and the 
magnetohydrodynamic Penrose process, we examine energy extraction from a 
black hole caused by magnetic reconnection in the ergosphere. The reconnection 
redistributes the angular momentum efficiently to yield the negative energy-at- 
infinity. We derive a condition for the process to operate in a simple situation, 
where the plasma is incompressible and the magnetic energy is converted com- 
pletely to the plasma kinetic energy locally. Astrophysical situations of magnetic 
reconnection around the black holes are also discussed. 

Subject headings: Black hole physics, magnetohydrodynamics: MHD, relativity, 
methods: analytical, galaxies: nuclei, gamma rays: bursts, plasmas 



1. Introduction 



Energy extraction from a rotating black hole interests us not only as engines of relativis- 
tic jets fr om active galactic nuclei (AGNs), micro-quasars (/iQSOs), and gamma-ray bursts 
(GRBs) (IMeier et al. 20011 ). but also as fundamentals of black hole physics. The horizon 
of the black hole is defined as the surface where no matter, energy, and information pass 
through outwardly. On the other hand, reducible energy from the rotating black hole is 
given by 



rot 



(1) 
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where M is the mass, a is the rotation parameter of the black hole, and c is the speed of 
light (IMisner et al. 1970l ). It corresponds to the rotational energy that can be extracted 
principally. 

Several kinds of distinct mechanisms have been proposed for the extraction of the rota- 



magnetohydrodynamic (MHD) Penrose process, and modified Hawking process ( 


Penrose 1969 




Press et al. 1972: 


Ruffini et al. 1975: 


Blandford et al. 1977: 


Hirotani et al. 1992: 


van Putten 2000) 



Here we mention three kinds of mechanism of the black hole energy extraction among them: 
the Penrose process, the Blandford- Znajek mechanism, and the MHD Penrose process, while 
the super-radiant scattering and the modified Hawking process may be related with the high 
energ y phenomena suc h as origins of gamma-ray bursts and the ultra-hig h energy cosmic 



rays (Ivan Putten 19991 : Ivan Putten 2000l : iPierre Auger collaboration 20071 ). The Penrose 



process involves production of particles with negative energy- at-infinity via strong fission 



or particle interaction in the ergosphere (jPenrose 19691 ). It needs drastic redistribution of 



angular momentum to produce the negative energy- at- infinity. Although this process clearly 
shows a possibility of energy extraction from the black hole, it is improbable as engines of 
astrophysical jets, because of poor collimation of particles and poor event rate. That is, 
the Penrose process accelerates the particles toward the equatorial plane, not toward the 
axis direction, and needs the azimuthal relativistic fission at not so wide region in the er- 
gosphere. Blandford & Znajek (1977) proposed a mechanism of energy extraction from a 
rotating black hole through force-free magnetic field. Their analytic steady-state solutions 
show the direct energy radiatio n from the horizon , which looks contradictory to the defini- 



tion of the black hole horizon (jPunsly et al. 1990l ). As we discuss in the next section, this 



mechanism also utilizes the negative energy-at-infinity. In this case, however, the negative 
energy-at-infinity is sustained not by the particle or matter, but by the electromagnetic 
field. Angular momentum of the electromagnetic field is redistributed by the magnetic ten- 
sion through almost mass-less plasma to produce the negative energy-at-infinity of the field. 
The magnetic tension may also redistribute the angular momentum of the plasma to yield 
the negative energy-a t-infinity of the plasma, when the plasma has non-zero mass density 



(IHirotani et al. 19921 ). It is called the MHD Penrose process. This energy extraction was 



confirmed by the n umerical simulations based on the general relativistic magnetohydro dy- 



namics (GRMHD) flKoide et al. 20021 : iKoide 2003h . 



It is noted that magnetic reconnection also redistributes angular momentum of the 
plasma to form the negative energy-at-infinity because it produces a pair of fast outflows 
with the opposite directions from the reconnection region. Then the rotational energy of the 
black hole can be extracted through the induced negative energy-at-infinity of the plasma. 
In the present paper, we derive a condition for the process to operate in a simple situation 
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for the incompressible plasma, where all the magnetic energy is converted to the plasma 
kinetic energy. 

In ^ we review the mechanisms of energy extraction from the rotating black hole in 
terms of the negative energy-at-infinity. In ^ we examine the operation condition of the en- 
ergy extraction from the black hole induced by magnetic reconnection in the ergosphere using 
a simple model. In ^ we discuss astrophysical situations where the magnetic reconnection 
happens in the ergosphere. 



2. Penrose process, Blandford-Znajek mechanism, and MHD Penrose process 

We use the Boyer-Lindquist coordinates {ct, r, 9, 0) to describe the space-time around a 
rotating black hole. The line element of the axisymmetric, stationary space-time around the 
black hole is written by 



(2) 



i=l 



where hi is the scale factor of the coordinate a;*, Ui is the angular velocity describing a frame- 
draggin g effect, a is the lap se function, and Pi = hiUJi/{ca.) is the shift vector. For the Kerr 
metric (IMisner et al. 1970l ). we have 



a 



AS 



, hi 



ho =v^, h-. 



— sin 6*, uji=uj2 = 0, uj3 = — -| — , (3) 



where A = - 2rgr + (arg)^, T, = r"^ + (arg)^ cos^ 6^, A = {r"^ + (arg)^}^ - A(arg)^ sin^ 61, 
and Tg = GM/c^ is the gravitational radius of the black hole. 

This metric has translational symmetry with respect to t and (p, so that we obtain the 
conservation law 



d 



(4) 



where g = Det{g^,y) = —{ahih2h^)'^ is the determinant of the metric tensor, T'^'^ is the 
energy-momentum tensor, and ^'^ is the Killing vector. When we adopt one component 
approximation of the plasma, we have 



T^"" = pg'''' + i)UW + F^F"'' - -g^'^^FP'^Fp^, 



(5) 



where p is the proper pressure, f) = Cint + p is the e nthalpy density, is the four- velocity, 
and Fp_jj is the electromagnetic field-strength tensor (IKoide et al. 19991 ). The thermal energy 
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density is given by cint = p/(T — 1)+ pc^ for adiabatic plasma, where F is the adiabatic index 
and p is the proper mass density. 

When we consider the Kilhng vectors C,'^ = (—1,0,0,0) and (0, 0, 0, 1), we get the 
energy and angular momentum conservation laws, 

de°° 1 ^ d 



h;-Xlur.{h.h,h,S^). (6) 



dt hih2hs ^ dx^ 

I - -MkE|r(M.M"), (7) 

i 

where e°° = —OLgy^T^^ is called "energy-at-infinity" density, which corresponds to the total 
energy density of the plasma and field, 5** = —cagyoT^" is the energy flux density, / = 
ag^vT^" jc is the angular momentum density, and M* = ahiT^^g^^ is the angular momentum 
flux density. 

When we introduce the local frame called the "zero angular momentum observer" 
(ZAMO) frame, we have di = adt, dx^ = hi{dx^ — ujidt). Because this is the local Minkowski 
space-time: ds^ = —{cdi)'^ + ^^^^{dx^y = rj^^dx^dx^ , the variables observed in the frame are 
intuitive. For example, the velocity -0*, the Lorentz factor 7 = f/" = [l — X]i=i(^Vc)^] ? 
and the four-velocity [/* = hiV^ — cajS^U^ {i = 1, 2, 3) have the relation, f/* = 7^*. Hereafter 
we denote the variables observed in the ZAMO frame with the hat, A. From equation 
we obtain 

e°° = ae + ^ cOihiP' = ae + cu^l, (8) 

i 

I = hP\ (9) 

Here e = P)7^ — p + [b^ + E'^/c^^ /2 is the total energy density and = {f)7^'0* + (^^ x 

By}/c'^ is the z-th component of the momentum density, where = and E"^ = 

o]^^ FjQ = cFq = cFio [i = 1, 2, 3) are the magnetic flux density and the electric field, 
respectively [e'^"^ is the Levi-Civita tensor). We normalize the field strength tensor F^ 



so that l-Bp/2 and |£'p/(2c^) give the magnetic and electric energy densities, respectively. 
For example, the magnetic field measured in the SI unit divided by the square root of the 
magnetic permeability po is B and the electric field measured in the SI unit times the 
square root of permittivity of vacuum eo is E. Equations ([8]) and IQ can be separated into 
the hydrodynamic and electromagnetic components: e°° = + eg^ and / = Zhyd + ^em, 
where 



e^d = "^hyd + '^ujih~-f^v' = achyd + ^^skyd, (10) 
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EM 



aCEM + ^i^i-^ ^ ^) . = "^EM + t^sW, (H) 



/hyd = /i3^7'^^', (12) 

Iem = h^(^ExBy^. (13) 

Here, Chyd = 1^7^ — P and cem = (-B^ + &/c^)/2 are the hydro dynamic and electromagnetic 
energy densities observed by the ZAMO frame, respectively. 



2.1. Penrose process 

When we consider a particle with rest mass m at rp{t) in the absence of electromagnetic 
field, i.e. p = ^S^{r — rp{t)), p = 0, and B = E = 0, then the energy- at- infinity and the 
angular momentum of the particle are given from equations ([H]) and as. 



E~ = / e'^dV = a-fmc^ + u;3L, (14) 



'n 

L = I IdV = hsm-^v^, (15) 
Jn 

where Q is the whole volume of the space. Both the energy- at- infinity and angular momentum 
of the particle conserve when it travels alone. 

Equations (fT^ and (fTSl) yield the energy-at-infinity of the particle as 

= a7mcM 1 + /33— ) . (16) 



£,3 
C 

-,3 



If Psv /c < —1, the energy-at-infinity of the particle becomes negative. This condition can 
be satisfied only in the ergosphere (/?3 > 1). Using the relation {1 — J^iiPiY} = "fi'oo, we 
have the well-known definition of the ergosphere: Qq^ > 0. 

When we consider the particle fission in the ergosphere, A — > B + C, the conservation 
laws of the energy-at-infinity and the angular momentum are 

= E^ + E^, (17) 
La = Lb + Lc, (18) 
= a-fimic^ + uj^Lj (/ = A,B,C). (19) 

If the fission is so strong that satisfies Lb = /^-3W^B7B'0B < —ch^rriB/ ■>/ (/^s)^ — 1, then we get 
the negative energy-at-infinity E^ < and E"^ < Eq . The particle C escapes to infinity 
and the particle B is swallowed by the black hole to reduce the black hole mass. Eventually, 
the rotational energy of the black hole is extracted. 
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2.2. Blandford-Znajek mechanism 



In the Blandford-Znajek mechanism, the energy propagates outwardly from the black 
hole horizon when the angular velocity of the black hole horizon is larger than that 
of the magnetic field lines Qp: fin > ^f- Although this statement looks inconsistent to 
the definiti on of the horizon, where any energy, matter, and information never transform 
outwardly (jPunsly et al. 19901 ) . it can be understood as the transportat ion of the n egative 
electromagnetic energy-at-infinity of magnetic fields into the black hole (IKoide 20031 ). Here 
we show the electromagnetic energy-at-infinity becomes negative when fin > ^f- 

From equations (ITTl) and (IT^ and the definition of cem and P\ the electromagnetic 
energy-at-infinity e^^j is written as 



-EM 



+ al3- [E X B 



(20) 



When we assume steady-state of the electromagnetic field in the force-free condition, the 
electric field observed by the ZAMO frame is given by 



E 



-^i^F 
a 



UJ3) e^x B 



(21) 



where e,^ is the unit vector parallel to the azimuthal coordinate, B = B — B^e^ is the 
poloidal magnetic field, and VLf is a constant along the magnetic fiux tube. If we use the 
velocity of the magnetic fiux tubes observed in the ZAMO frame, Vf = {h2,/a)(yLF — i^a)^^, 
equation fl^ can be written as a intuitive equation, E = —Vf x B . Inserting equations 
( !2T|) into (!20|) . we have the electromagnetic energy-at-infinity. 



-EM 



a- 



P\2 



ca 



(/93 



(22) 



We consider the electromagnetic energy-at-infinity very near the horizon. We adopt the 
boundary condition of the electromagnetic field at the horizon. 



B4, _ Vpjr th) 
~ c 



(23) 



where th is the radius of the horizon. The condition (!23l) is intuitive when we assume the 
Alfven velocity is the light speed c in the force-free magnetic field. It is also identical to 
the condition used by Blandford & Znajek (1977), while the original expression is rather 
complex. After some manipulations with equations fl22l) and fj23l) . we obtain at the horizon 

2 

(5P)21]f(1^f-^3). (24) 
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Here we have used a — * very near the horizon. The electromagnetic energy-at-infinity 
is negative only when < l^p < c^s, which is identical to the switch-on condition of the 
Blandford-Znajek mechanism. The coincidence indicates that the Blandford-Znajek mech- 
anism utilizes the negative electromagnetic energy-at-infinity. This conclusion is generally 
applicable to the force-free field of any spinning black hole as far as equation (123|1 is valid, 
while the original analytical model of Blandford & Znajek (1977) is restricted to a slowly 
spinning black hole. In the force-free condition, the Alfven velocity becomes the speed of 
light, and then the Alfven surface is located at the horizon. Then the region of the negative 
energy-at-infinity can be connected with the region outside of the ergosphere through the 
magnetic field causally. 



2.3. MHD Penrose process 

The MHD Penrose process is the mechanism of energy extraction from a black hole 



through the negative energy-at-infinity of plasma induced by the magnetic tension (IHirotani et al. 19921 ) 
It has been confirmed by Koide et al. (2002) and Koide (2003) using GRMHD numerical sim- 
ulations. It is defined as the energy extraction mechanism with the negative energy-at-infinity 
of plasma, which is induced by the magnetic tension, while in the Blandford-Znajek mech- 
anism, the negative electromagnetic energy-at-infinity plays a important role. For a rapidly 
rotating black hole (a = 0.99995), the magnetic fiux tubes of the strong magnetic field which 
cross the ergosphere are twisted due to the frame- dragging effect. The angular momentum of 
plasma in the ergosphere is opposite to that of the rotating black hole and its magnitude is 
large to make negative in equation (fTOj) . The twist of the magnetic flux tubes propagates 
outwardly. The Poynting flux indicates that the electromagnetic energy is radiated from the 
ergosphere (see flgure 4 of Koide et al. 2003). At the foot point of the energy radiation from 
the ergosphere, the hydrodynamic energy-at-inflnity ej^^ = a7ehyd + ^^s^hyd decreases rapidly 
and becomes negative quickly. The negative energy-at-inflnity is mainly composed of that 
of the plasma. To realize the negative energy-at-inflnity of the plasma, redistribution of the 
angular momentum of the plasma, Zhyd = h^'^'^'^v^, is demanded. The angular momentum of 
the plasma is mainly redistributed by the magnetic tension. 



3. Magnetic reconnection in ergosphere 

We investigate energy extraction through negative energy-at-inflnity induced by mag- 
netic reconnection in an ergosphere around a Kerr black hole. For simplicity, we consider the 
magnetic reconnection in the bulk plasma rotating around the black hole circularly at the 
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equatorial plane (Fig. [T]). To sustain the circular orbit, the plasma rotates with the Kepler 
velocity fx or it is supported by external force, like magnetic force. Here the Kepler velocity 
is given by 



The plus (minus) sign corresponds to the co-rotating (counter-rotating) circular orbit case. 
Throughout this paper, we use the co-rotating Kepler velocity. We assume the initial anti- 
parallel magnetic field directs to the azimuthal direction in the bulk plasma, and the pair 
plasma outflows moving toward the opposite directions each other caused by the magnetic 
reconnection are ejected in the azimuthal direction. It is also assumed that the plasma ac- 
celeration through the magnetic reconnection is localized in the very small region compared 
to the size of the black hole ergosphere, and the magnetic field outside of the plasma acceler- 
ation region is so weak that the plasma flow accelerated by the magnetic reconnection is not 
influenced by the large-scale magnetic field around the black hole. If one of the pair plasma 
flows in the opposite direction of the black hole rotation has negative energy-at-infinity and 
the other in the same direction of the black hole rotation has the energy-at-infinity higher 
than the rest mass energy (including the thermal energy) (Fig. [T]), the black hole rotational 
energy will be extracted just like the Penrose process. 

We have to investigate two conditions, i.e. the condition for the formation of the 
negative energy-at-infinity and the condition for escaping to infinity. Before we move on to 
the conditions, we here mention about the elementary process of the relativistic magnetic 
reconnection in the locally uniform, small-scale plasma rotating circularly around the black 
hole. To investigate the magnetic reconnection in the small scale, we introduce the local 
rest frame {ct' of the bulk plasma which rotates with the azimuthal velocity 

= cf3o at the circular orbit on the equatorial plane r = < r^, 6 = tt /2. We set the frame 
[ct', x^', x"^', x^') so that the direction of x^' coordinate is parallel to the radial direction x^ = r 
and the direction of x^' is parallel to the azimuthal direction x^ = (p (Fig. [I]). Hereafter we 
note the variables observed by the rest frame of the bulk plasma with the prime, "/". First, 
we consider the magnetic reconnection in the rest frame locally and neglect the tidal force 
and Coriolis' force for simplicity. We regard that the rest frame rotating with the Kepler 
velocity is in a gravity-free state and thus we can consider the magnetic reconnection in 
the framework of special relativistic MHD. The initial condition is set to that of the Harris 
model where the anti-parallel magnetic field and the plasma are in equilibrium: 



cA -^^pr^jr — ary r"^ 




(25) 




B 



3/ 



Sotanh(x^7(5), B 



B 



21 



(26) 



V = 




(27) 
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V 



P 



Po, 







,2/ 







V' 



(28) 
(29) 



where Bq is the typical magnetic field strength and 26 is the thickness of the current layer 
(see Fig. [2]). The electric resistivity is assumed to be zero except for the narrow reconnection 
region. We assume the system is symmetric with respect to the a;^'-direction, and resistivity 
is given by 



where t]q is a positive constant and f{x^',x^') is the profile of the resistivity, which is finite 
around x^' = 0, x^' = but zero outside the reconnection region. The magnetic fiux tubes 
reconnected at the resistive region accelerate the plasma through the magnetic tension, and 
the magnetic energy in the fiux tubes is converted to the kinetic energy of the plasma. As 
shown in Fig. [3], the magnetic fiux tube and the plasma run away and fresh tubes and the 
plasma are supplied to the reconnection region from outside of the current layer successively. 

The outfiow velocity of the accelerated plasma through the magnetic reconnection, f^^^., 
is estimated by the velocity v'^^^, where whole magnetic energy is converted to the kinetic 
energy. If the magnetic energy is completely converted to the kinetic energy of the plasma 
particles which are initially at rest, the magnetic energy per plasma particle is So/(2no), 
where no is the plasma particle number density. Using the approximation that the plasma 
element is treated as incompressible gas covered by very thin, light, adiabatic skin with the 
thermal energy U and the enthalpy H, using equation (lASp in Appendix A, we can write 
the energy conservation equation with respect to the particle as. 



where 7^3^^ is the Lorentz factor of the plasma particle after complete release of the magnetic 
field energy observed in the bulk plasma rest frame. We get the maximum Lorentz factor 
y^^yr from equation fl3T|) as. 



where u\ = B^/ho, zu = po/ho, AD = [u\ + 4:{l — w)]'^ + 16zu, po = (T — l)UnQ is the pressure, 
and ho = uqH is the enthalpy density of the plasma. Obviously, the terminal velocity of the 
plasma outfiow through the magnetic reconnection is smaller than the maximum velocity 
v'^^^, since all the magnetic energy is not always converted to the kinetic energy because of 
the Joule heating in the reconnection region and finite length acceleration. Figure H] shows 
the four-velocity u'^^^ = 7out'^out of the plasma outfiow caused by the magnetic reconnection 




(30) 




(31) 




(32) 



- 10 - 



against the plasma beta /3p = 2pq/Bq in the case of F = 4/3 and Bq 
denotes the predicted values from equation fl32l) . 



Pqc^- The solid line 



where 



7n 



Ua 



v' 1 

max 



1 + 



po(_ 

^0 



+ 



r 



2(r - 1 



-1/2 



+ 2(zu-r 



u. 



2 



+ 



r 



2(r- 1 



(33) 



(34) 



The full squares are from the numerical calculations (Watanabe et al. 2006). It is found 
that our values of wJnax good agreement with the result of numerical simulations. In 

the high plasma beta region, the values of u'^^^ from equation fl33l) are slightly smaller than 
those of the numerical result. This discrepancy, i.e. M^^ax < ^out^ should come from the 
release of the thermal energy to the plasma kinetic energy. Therefore, we take the maximum 
velocity v'^^^^ as the plasma outflow velocity induced by the magnetic reconnection hereafter. 

Now we return to consider the conditions of the negative energy-at-infinity formation 
and escaping to infinity of the pair outflows caused by the magnetic reconnection. From 
equations (1A6I) . the hydrodynamic energy-at-infinity per enthalpy of the plasma ejected 
through the magnetic reconnection into the ±x^' (azimuthal) direction with the Lorentz 
factor 7^^^^ is given by 



e± (MA,/3o,tJ7,Q;,^3) = 



-hyd,± 



"70 



:i + /?3/9o)7;ax ± 00 + /^S) V(7;nax)' " 1 



7maxT/3oy(7 



' )2 - 1 

max/ 



(7U)^ + m 



-w 



(35) 



where 7^3^^ is given by equation ( l32i) . The plus and minus signs of the subscript in express 



the cases with the plasma velocity v 



3/ 



and V 



3/ 



~^max' respectively. If < 0, 



the energy-at-infinity of the plasma becomes negative. Here we neglect the contribution 
from the electromagnetic field, because when the plasma velocity becomes fmax; most of 
the magnetic energy is converted to the plasma kinetic energy so that the electromagnetic 
energy-at-infinity becomes negligible. Then the total energy-at-infinity of the plasma and the 
electromagnetic field becomes negative. Here e'^ decreases monotonically when increases, 
and it is positive when ua = and negative when ua is large enough in the ergosphere. Then, 
we can define the one- valued function f/^(/3o, w, a, j3z) which satisfies ef°(f/^, /3o, w, a, (3^) = 
0. The condition for the formation of the negative energy-at-infinity through the magnetic 
reconnection is given by > U^. 



The condition of escaping to infinity of the plasma particle accelerated by the magnetic 
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reconnection is given by 

Ae- ^er-(l- ^t^) > 0. (36) 

Here we assume the magnetic field is so weak far from the reconnection region that we 
can neglect the interaction between the large-scale magnetic field and the outflow from the 
reconnection region. The monotonic function Ae^^ with respect to ua is negative when 
ua = and positive when ua is large enough. Then we can also define the one-valued 
critical function Uj^0o,w,a, P^j such that Ae'^{Up^, Po,uj,a, Ps) = 0. The condition of the 
escape to infinity of the plasma flow is written by ua > U^. 

We apply the conditions ua > and ua > to the plasma with the rotation bulk 
velocity -0^ = c/3o = ^k- Figure [5] shows and against r/rs (rs = 2rg) for the fixed black 
hole rotation parameter a = 0.995 and the pressure parameter w = 0, 0.2, 0.5. The lines of 
the critical values of Alfven four-velocities and are drawn as the contours of e^^ = 
and AeJ^ = 0, respectively. The vertical thin dashed line indicates the inner most stable 
orbit radius of a single particle. The vertical thin dot-dashed line shows the point of -Ok = c. 
Then, the Kepler motion is unstable between the dot-dashed line and the dashed line. In 
the left region of the dot-dashed line, there is no circular orbit anymore. The vertical thin 
solid line indicates the horizon of the black hole. The upper regions of the thick lines e!^ = 
show the condition of the formation of the negative energy-at-infinity through the magnetic 
reconnection, ua = -Bo/v^ > U^. The different styles of lines indicate the different pressure 
cases (solid line: zu = 0, dashed line: vj = 0.2, dot-dashed line: zu = 0.5). The regions above 
the thick lines Ae^^ = show the condition on escaping of the plasma accelerated by the 
magnetic reconnection, ua = Bo/y/ho > U^. The difference of the line styles denotes the 
same as that of the upper lines. It is shown that the condition of the energy extraction 
from the black hole through the magnetic reconnection is determined from the condition of 
the formation of the negative energy-at-infinity, ua > in the case of the rotating plasma 
with the Kepler velocity. In the zero pressure case (ci7 = 0), the easiest condition is found 
at r = 0.61rs, ua > = 0.86. This means the relativistic reconnection is required for the 
energy extraction from the black hole. In the finite pressure case {w = 0.2), the condition is 
relatively relaxed around the outer region of the ergosphere, while the condition around the 
inner ergosphere is severer. However, the difference is small between these cases. This case 
also requires the relativistic magnetic reconnection to extract the black hole energy. The 
condition of the case w = 0.5 is almost similar to the tendency of the previous two cases 
(ro = 0,0.2). 

Figure [6] shows the critical Alfven four- velocity of the energy extraction from the 
black hole with the rotation parameter a = 0.9. In this case, the condition of the energy 
extraction from the black hole becomes severer compared to the cases with the larger rotation 
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parameters a = 0.995. 

It is noted that in the case of the black hole rotation parameter a < 1/^/2, there is 
no circular orbit inside of the ergosphere. In such case, we can not consider the energy 
extraction through the magnetic reconnection in the circularly rotating plasma around the 
black hole except for the case with the support by magnetic field. 

When we consider the slower rotating plasma case c(3o < vk, which may be an artificial 
assumption, the conditions of the negative energy-at-infinity and the escape plasma become 
comparable and relaxed as a whole. Figure [7] shows the critical Alfven four-velocity of the 
case c/3o = 0.6vk, a = 0.995. At r = 0.6rs, ~ f/^ ~ 0.5 for both the tu = and 
w = 0.2 cases. This means the sub-relativistic magnetic reconnection can extract the black 
hole energy. 

4. Discussion 

In the previous section, we showed the possibility of the energy extraction from the black 
hole through the magnetic reconnection in the ergosphere. In this mechanism, the magnetic 
tension plays a significant role to cause the plasma flow with the negative energy-at-infinity, 
like the MHD Penrose process. If we consider quick magnetic reconnection, the mechanism 
by the magnetic reconnection is more effective than the MHD Penrose process, because the 
fast plasma flow can be induced, so that all magnetic energy can be converted to the kinetic 
energy of the plasma flow through the magnetic reconnection. 

As the magnetic reconnection mechanism, we utilized a rather simple model with an 
artificial resistivity at the reconnection region and the local approximation around the rotat- 
ing black hole. We further assumed that the outflow caused by the magnetic reconnection 
is parallel to the azimuthal direction. In general, the outflow is oblique to the azimuthal 
direction. In the oblique case with the angle x between the outflow and azimuthal direc- 
tions, the condition of the energy extraction through the magnetic reconnection is given by 
Ma > U^{ua,'!^,(x, Pscosx) and ua > f^A (""a, "W, «, /^s cosx), where of equation (135|) is 
replaced by /^s cos x- This is a severer condition compared to the parallel case. 

Let us briefly estimate the critical magnetic field required for the energy extraction from 
the black hole through the magnetic reconnection with respect to the AGNs, /xQSOs, and 
GRBs. Here the critical magnetic field in the SI unit, -Bcrit, is given by i?crit = V/^o^o ~ 
^fioPoC, where Hq and po are the typical enthalpy and mass density around the objects, 
respectively. Here, we assume the pressure is not larger than poc^ and neglect it to get rough 
estimation. To estimate the typical mass density po of the plasma around the central black 
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hole of these objects, we use 



Po 



5 X 10^ 



M 



M 



2r 



-3/2 



gem , 



(37) 



where M is the accretion rate and M is the black hole mass (IRees 19841 ) . For the AGN in the 
l arge elliptical galaxy M8 7, whe n we assume M = 10~^Mq yr~^, r = r^, and M = 3 x IO^Mq 
( IRevnolds et al. 19961 : IHo 19991 ). equation (l37Il yields po — 2 x lO^^'^gcm"^. Then, we get 
Bcrit — 500 G. This magnetic field is probable around a black hole of an AGN, thus the 
magnetic extraction of black hole energy is possible. 



In a collapsar model with M = O.IMqS"^ and M = 3Mq flMacFadven et al. 19991 ) 



when we apply equation fl37j) . we get po — 6 x lO^gcm ^ at r = rs as mass density around a 
black hole in a GRB progenitor. The critical magnetic field is then Brr,t 8 x 10^^ G . The 



magnetic field of the progenitors is estimated to be lO^^G to lO^^G ( Ivan Putten 19991 ) and 
then extraction of the black hole energy through the magnetic reconnection is marginally 
probable in a core of a collapsar. 

With respect toyxOSO, GRS1915+1 05 has a mass accretion rate of M = 7x 10~'^Mm yr~ ^ 
jMirabel et al. 19941 : IPender et al. 20041 ) with a mass of M = 14M0 JGreiner et al. 200ll ). 



Then equation fl37|) yields po ~ 6 x 10~^gcm~^. The critical magnetic field is estimated as 
Bcrit — 8 X 10*^ G. This magnetic field is too strong as the field around a black hole in pQSOs. 
Thus, the energy extraction from a black hole in pQSO through the magnetic reconnection 
may not be prospective. 

We discuss the possibility of formation of anti-parallel magnetic field with a current 
sheet where strong magnetic reconnection in the ergosphere is caused. First, let us con- 
sider uniform magnetic field around a rotating black hole as the initial condition. In this 
magnetic configuration, one may think that the magnetic reconnection scarcely happens. 
However, this is caused naturally by the gravitation and the frame-dragging effect of the 
rapidly rotating black hole. Under this situa t ion, GRMHD simulations we re carried out 
with zero electric resistivity (iKoide et al. 20021 : iKoide 20031 : iKomissarov 20041 ). The numer- 
ical simulations showed that the magnetic fiux tubes across the ergosphere are twisted by 
the frame-dragging effect of the rotating black hole, and the plasma falling into the black 
hole makes the magnetic field radial around the ergosphere. Here it is noted that the mag- 
netic field line twisted by the frame-dragging effect makes the angular momentum of the 
plasma around the equatorial plane and the ergosphere negative (/ < 0), and the plasma 
with the negative angular momentum falls into the black hole more rapidly. The attractive 
force toward the black hole comes from the shear of the frame-dragging effect directly (e.g. 
see the term with aji in equation (56) of Koide (2003)). Beside the equatorial plane in the 
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ergosphere, the magnetic field becomes anti-parallel. The magnetic flux tubes are twisted 
strongly enough, and then the strong anti-parallel open magnetic fleld is formed almost along 
the azimuthal direction (Fig. [H^). In this way, the magnetic reconnection happens and the 
energy of the rotating black hole is extracted through the magnetic reconnection, even in 
the case of the initially uniform magnetic fleld. When the magnetic reconnection happens 
around the anti-parallel magnetic fleld, the outward flow from the reconnection region will 
be ejected toward inflnity along the open magnetic fleld lines. This outflow will be bent and 
pinched by the magnetic fleld and may become a jet. 

Next, as the initial condition, we assume closed mag netic flux tubes which are believed to 



be formed in the accretion disks around the black holes (Ivan Putten 19991 : iKoide et al. 2006 



McKinney 20061 ) (Fig. [Hb). When a single closed magnetic flux tube is tied to an edge of a 
rotating quasi-stationary disk and a bulk part of the disk, the plasma at the edge loses the 
angular momentum and falls into the black hole, while the bulk plasma tied to the magnetic 
flux tube increases the angular momentum and shifts outwardly. The plasma at the disk 
edge falls spirally due to the frame-dragging effect, and the magnetic flux tubes dragged by 
the plasma are elongated spirally. If the twist of the magnetic flux tube is strong enough, 
anti-parallel closed magnetic fleld is formed almost along the azimuthal direction. In this 
magnetic conflguration, energy may be extracted from the black hole through the magnetic 
reconnection in the ergosphere. The outflow from the reconnection region will elongate the 
closed magnetic fleld lines, while this decelerates the outflow. If the magnetic reconnection is 
caused in the elongated magnetic fleld lines, the plasmoid is formed and is ejected to inflnity. 

As shown in the above two cases of the open and closed magnetic fleld, the large- 
scale dynamics of the outflow through the magnetic reconnection depends on the large-scale 
magnetic field configuration. Here we note that closed magnetic fiux tubes across an accretion 
disk and an ergosphere around a rapidly rotating black hole is unst able and expands vertically 



to form open magnetic field (IKoide et al. 20061 : iMcKinney 20061 ). Then, around a rapidly 



rotating black hole, open magnetic field may be probable compared to closed field. Anyway, 
these phenomena should be investigated by numerical simulations of the full GRMHD with 
non-zero electric resistivity (resistive GRMHD). 

To be more exact, in both cases of the open and closed topologies of the magnetic field, 
the outfiow caused by the magnetic reconnection is infiuenced by the large-scale magnetic 
field, and it is not determined only from the local approximation which we used here. In 
both cases, it is noted that the plasma with the negative energy-at-infinity is farther from 
the horizon than the plasma accelerated by the magnetic reconnection. If we assume the 
azimuthal symmetry of the initial condition, interchange instability should be caused so 
that the plasma with negative energy-at-infinity falls into the black hole and the plasma 
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with additional energy- at- infinity runs away to infinity. With respect to the interchange 
instabihty, we consider essentially hydrodynamic mode where initially super-Keplerian inner 
part of the disk supports sub-Keplerian outer part against the black hole gravity. The 
closed magnetic flux tube is formed across the plasma with additional energy-at-infinity. 
The plasma at the inner edge of the magnetic loop falls into the black hole because of the 
deceleration by the magnetic tension, while the plasma at the outer edge of the magnetic 
loop is accelerated and escapes to infinity. Then the magnetic flux tube is elongated between 
the escaping and falling plasmas. In such magnetic flux tube, the anti-parallel magnetic field 
with strong current sheet may be formed and the magnetic reconnection may be caused 
once again. On the other hand, the plasma with the negative energy-at-infinity through 
the magnetic reconnection falls into the black hole and the magnetic flux tube tied to the 
plasma is also elongated by the frame-dragging effect. The anti-parallel magnetic field in 
the magnetic flux tube will also form and the magnetic reconnection is caused repeatedly. 
Above discussion shows that the magnetic reconnection can be caused intermittently in the 
ergosphere. To investigate these phenomena, the numerical simulations of resistive GRMHD 
are also demanded. 

The resistive GRMHD should solve the problems with respect to the energy extraction 
through the magnetic reconnection. For example, using resistive GRMHD simulations, we 
can take into account of the situation that the initial plasma falls into the black hole with sub- 
Keplerian velocity of the plasma rotation, which is neglected in this paper. Unfortunately, 
no simulation with resistive GRMHD has been performed until now, while recently, ideal 
GRMHD numerical simulations, where the electric resistivity is zero, have come spread 
(Koide et al. 1998, 1999, 2000, 2002, 2006; Koide 2003, 2004; McKinney 2006; Punsly 2006; 
Komissarov et al. 2007, and references therein). They confirmed important, interesting 
magnetic phenomena around the rotatin g black holes, such as magnetically-induced energy 



extraction from the rotating black hole (IKoide et al. 20021: iKoide 20031 : iKomissarov 20041 ) 



and formation of magnetically-driven relativistic jets (IMcKinney 20061 ). To confirm that 
the plasma with additional energy by the magnetic reconnection escapes to infinity and the 
plasma with the negative energy-at-infinity falls into the black hole, we may also use the ideal 
GRMHD simulations. In such calculations, we can trace the plasma trajectories after the 
magnetic reconnection stops. The magnetic configuration of the post stage of the magnetic 
reconnection is used as an initial condition. 

On the other hand, in spite of the restriction of th e ideal GRMHD, many magnetic is- 



lands are seen in the last stages of long-term calculations (IMcKinney et al. 20041 : iKoide et al. 2006 



McKinney 20061 ). These magnetic islands, of course, are artificial appearance. However, 
these numerical results indicate that the magnetic configuration where the magnetic recon- 
nection occurs is relatively easily formed around the black hole. Recent X-ray observations 
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of the solar corona confirmed that the magnetic reconnection takes places frequently in the 
active region of the corona and causes drastic phenomena, like solar fiares. The observations 
and the recent theories of solar and stellar fiares indicate that the magnetic reconnection 
is common in the ast rophysical plasmas around t he Sun, stars, and the compact objects 
including black holes (IMasuda 1994 IShibata 19971 ). The energy extraction from the black 
hole through the magnetic reconnection is one of the phenomena of magnetic reconnection 
around the black hole. More drastic phenomena related with the magnetic reconnection may 
exist. To investigate these phenomena, resistive GRMHD numerical calculations will play a 
crucial role. 
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A. Relativistic adiabatic incompressible ball approach 

To take account of inertia effect of pressure into plasma, we use an approximation of 
incompressible fluid, which is consisted of small, separated, constant volume elements. The 
element is covered by thin, light, adiabatic, closed skin and its volume is constant, like a ball 
used for soft tennis. We call this method "relativistic adiabatic incompressible ball (RAIB) 
approach". Here we assumed gas pressure does not work to the plasma and influences the 
plasma dynamics only through an inertia effect. Let us consider the fluid in one ball with 
the mass m. When the ball locates at r = r(t), the mass density of the gas is 

TT) 

p{r,t) = —-6^r-r{t)), (Al) 

lit) 

where 7(t) is the Lorentz factor of the ball at time t and 6'^{r) is the Dirac's (5-function in 
three-dimensional space. Because the gas in the ball is assumed to be incompressible and 
adiabatic and then its temperature is constant, the pressure should be proportional to the 
mass density, p{r,t) oc p{r,t). On the other hand, the thermal energy in the ball 

U = l^ ^-^l{t)dV (A2) 
should be constant. Then we found 

(r-i)f/,3, 



p{r,t) = ^^^6\r-r{t)). (A3) 
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Using equations ffTOl) . flAll) . and flA3|l . we found the energy-at-infinity of the gas in the ball 

as, 



dV 



a 



-P 



dV 



a 



(A4) 



where (3^ = v^/c and H = rac^ + TU. In the special relativistic case, equation flA4p yields 
the total energy of the ball as. 



tot 



7iJ 



-U. 



7 



(A5) 



Next we derive the energy-at-infinity of the incompressible ball rotating circularly around 
the black hole. We assume that the bulk plasma rotates circularly with the velocity -O^ = c/3o, 
and then the relative three- velocity between the rest frame of the bulk plasma and the ZAMO 
frame is v'^ = c/3o, = = 0. The line elements of the ZAMO frame {ct ^) and 

the bulk plasma rest frame (ct', x^', x^', x^') are related by the Lorentz transformation. Using 
equations flA4p and the Lorentz transformation, the energy-at-infinity of the incompressible 
ball with four- velocity observed by the bulk plasma rest frame (7', U'^, U'^, U'^) is given by 



:,2/53 



7" + lif^l 



(A6) 



It is noted that the four-velocity of the plasma is related with the Lorentz factor by f/'^ = 
±v/y2"^ in the case of U'^ = U'"^ = 0. This formula of the energy-at-infinity of the 
incompressible ball ( ]A6I) is applicable to that of one particle of the plasma effectively. 
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Fig. 1. — Schematic picture of the extraction of the black hole energy through the magnetic 
reconnection in the bulk plasma which rotates circularly on the equatorial plane of the 
rotating black hole. The phenomena are indicated in the Boyer-Lindquist coordinates. The 
coordinates O' — x^'x^'x^' in the inserted box are the rest frame of the bulk plasma. 
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Fig. 2. — Initial magnetic configuration around the reconnection region in the local rest 
frame of the bulk plasma. The area in the dark gray ellipse at the origin corresponds to the 
magnetic reconnection region and the gray layer along the x^' shows the current sheet. 




Fig. 3. — Energy convergence from the magnetic energy to the kinetic energy through the 
magnetic reconnection in the rest frame of the bulk plasma. 
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Fig. 4. — Comparison between the numerical result of the four-velocity of the plasma outflow 
Uo^t through the magnetic reconnection (full squares; Watanabe et al. 2006) and the simple 
expression (|33|) (solid hne) as a function of the plasma beta, /?p. 



- 23 - 




Fig. 5. — The critical Alfven four-velocity for the energy extraction from the rotating 
black hole with the rotation parameter a = 0.995 induced by the magnetic reconnection. The 
base plasma rotates around the black hole with the Kepler velocity c/3o = vk- The thick solid, 
dashed, and dot-dashed lines correspond to the cases of Ci7 = 0, 0.2, 0.5, respectively. The 
lines with the sign Ae'^ = {e'^ = 0) show the critical Alfven four-velocity of the condition 
on escaping to infinity of the plasma accelerated by the magnetic reconnection (the negative 
energy-at-infinity plasma formation). The vertical thin solid line at th = 0.550rs shows the 
black hole horizon. The vertical thin dashed line at r^s = 0.671rs indicates the radius of 
the marginal stable orbit. The vertical thin dot-dashed line at tl = 0.559rs shows the point 
where the Kepler velocity is the light velocity c. Between the vertical thin solid line and the 
vertical thin dot-dashed line, there is no circular orbit of a particle. Generally speaking, 
is infinite at r = rs, while is finite at this point. 
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Fig. 6. — Similar to Fig. O but for the black hole with the rotation parameter a = 0.9. The 
critical radii are rn = 0.718rs, tl = 0.779rs, and r^s > '"s- 
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Fig. 8. — Astrophysical magnetic configuration of tfie magnetic reconnection in tlie black 
liole ergospliere. (a) Anti-parallel magnetic field caused from the initial uniform magnetic 
field around the rapidly rotating black hole. The solid/dashed line shows the anti-parallel 
magnetic field line in front of /behind the equatorial plane from the reader. In this mag- 
netic field configuration, the current sheet locates at the equatorial plane, (b) Anti-parallel 
magnetic field formed by the closed magnetic flux tube tied to the disk around the rapidly 
rotating black hole. The inner part of the magnetic flux tube falls into the black hole to 
elongate the flux tube, and the anti-parallel magnetic field configuration is formed. 



